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Abstract 

A class of models is introduced describing the evolution of population species whose 
carrying capacities are functionals of these populations. The functional dependence of 
the carrying capacities reflects the fact that the correlations between populations can 
be realized not merely through direct interactions, as in the usual predator-prey Lotka- 
Volterra model, but also through the influence of species on the carrying capacities of 
each other. This includes the self-influence of each kind of species on its own carrying 
capacity with delays. Several examples of such evolution equations with functional 
carrying capacities are analyzed. The emphasis is given on the conditions under which 
the solutions to the equations display extreme events, such as finite-time death and 
finite-time singularity. Any destructive action of populations, whether on their own 
carrying capacity or on the carrying capacities of co-existing species, can lead to the 
instability of the whole population that is revealed in the form of the appearance of 
extreme events, finite-time extinctions or booms followed by crashes. 
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1 Brief survey of population models 



Evolution equations, describing population dynamics, are widely employed in various branches 
of biology, ecology, and sociology. The main forms of such equations are given by the vari- 
ants of the predator-prey Lotka-Volterra models. In this paper, we introduce a novel class 
of models whose principal feature, making them different from other models, is the func- 
tional dependence of the population carrying capacities on the population species. This 
general class of models allows for different particular realizations characterizing specific cor- 
relations between coexisting species. The functional dependence of the carrying capacities 
describes the mutual influence of species on the carrying capacities of each other, including 
the self-influence of each kind of species on its own capacity. Such a dependence is, both 
mathematically and biologically, principally different from the direct interactions typical of 
the predator-prey models. Before formulating the general approach, we give in this section a 
brief survey of the main known models of population dynamics. This will allow us to stress 
the basic difference of our approach from other models used for describing the population 
dynamics in biology, ecology, and sociology. 

(i) Predator-prey Lotka-Volterra model 

The first model, describing interacting species, one of which is a predator with population 
Ni, and another is a prey with population N 2 , has been the Lotka-Volterra [TJ[2] model 

^ = - Tl JV! + A 12 N 2 Nt , = l2 N 2 - A 21 N 1 N 2 , (1) 

at at 

where all coefficients are positive numbers. It is easy to show that the solutions to these 
equations are bound oscillating functions of time. 

(ii) Predator-prey Kolmogorov model 

The Lotka-Volterra model is a particular case of the predator-prey Kolmogorov model 
[3JH] that has the general form 

^ = fi(N 1 ,N 2 )N 1 , ^ = f 2 ( Nl ,N 2 )N 2 , (2) 
at at 

under the conditions 

> o «A < o 

8N 2 ' dNj 

This model is too general and requires specifications for describing concrete cases. 

(iii) Generalized predator-prey Lotka-Volterra model 

Generalizing the Lotka-Volterra model (QQ) to multiple species yields the equations 



dNi 
~dT 




(3) 



where all coefficients are real numbers [5]. The signs of the parameters can be different. 
When all 74 's are positive, while all A^s are negative, one gets the competitive Lotka- 
Volterra equations whose behavior has been analyzed in detail in Refs. [&H9]. 
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(iv) Replicator equations 

These equations have the form 



{fi -f)N t , (4) 



dt 

where fi is the species fitnesses and 

7=Yf* N * 

i 

is the average fitness characterizing the whole society [5]. The species populations are usually 
assumed to be defined on a simplex, being normalized to a constant representing the total 
fixed population 

^2 N t = N = const . 

i 

The n- dimensional replicator model is equivalent to the n — 1-dimensional Lotka-Volterra 
model ([3]), to which it can be reduced by a change of variables [51 [TU]. 

(v) Jacob-Monod equations 

The equations describe not the coexisting species but rather a single type of species of 
population iV\, like bacteria, which are fed on the nutrient of amount N 2 . The nutrient 
plays the role of the prey that is getting depleted being consumed by the feeders and, at 
the same time, being supplied into the system from outside according to the supply function 
f(N 2 ) = aN 2 /((3 + N 2 ). The equations read as 

^ = f(N 2 )m , ^ = -7/(JV 2 )JVi, (5) 

with all parameters being positive [UJ. As time increases, t — > oo, the nutrient becomes 
depleted, N 2 — > 0, and the bacteria population reaches the stationary value Ni = Ni(0) + 
iV 2 (0)/ 7 . 

The Holling equation of second kind [12] takes into account that predators, in order to 
consume prey, need to search for it, chase, kill, eat, and digest. This is why predators attack 
not all preys but a limited number of them, which saturates to a constant when the prey 
density increases. Mathematically, the Holling equation is analogous to the Jacob-Monod 
model. 

(vi) Verhulst logistic equation 

The well known logistic equation 

where all parameters are positive, was suggested by Verhulst [13]. The constant K is the 
carrying capacity. The solution to this equation is the sigmoid function 

N(t) " 



K + N (eT* - 1) 
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in which iVo = N(0). 

(vii) Hutchinson delayed logistic equation 

If one interprets the term inside the brackets in Eq. 



then, as Hutchinson argued 



) as an effective reproductive rate, 
it could be delayed in time, which results in the equation 



dN(t) 
dt 



7 JV(t) 



N(t 



T 



K 



(7) 



in which K is a fixed carrying capacity The solution to this equation gives additional 
oscillations superimposed on the logistic curve. 

(viii) Generalized delayed logistic equations 

There are many variants generalizing the delayed logistic equation (jTJ), which can be 
found in Refs. [ToT - fT?] . For example, the multiple-delayed equation 



dN(t) 
dt 



E 



N(t - TV) 



(8) 



where the carrying capacities Kj are positive constants. All such equations are usually 
applied to single-species systems of population N. The multiple carrying capacities Kj in 
Eq. (JSJ) correspond to different processes of the same single species. The logistic equations of 
the above type, whether with delays or without them, do not describe the possible coexistence 
of several species. When such equations are generalized to the case of several species, one 
comes back to the generalized Lotka-Volterra predator-prey model (]3]). 

(ix) Peschel-Mende hyperlogistic equation 

In order to take into account the accelerated growth of population, occurring, for instance, 
for the human world population, Peschel and Mende [18] extended the standard logistic 
equation by introducing two additional positive powers m and n, getting the equation 



N 

1 -K 



(9) 



The solution to this equation could be reasonably well fitted to the world population dy- 
namics. The form of the solution is a slightly modified sigmoid curve, with the population 
never surmounting the carrying capacity K. For m = 1 and n — 1, the Verhulst logistic 
equation (jSJ) is recovered. 

(x) Hyperlogistic time-delay equations 

The straightforward extension of the Peschel-Mende hyperlogistic equation is the hyper- 
logistic time-delay equation 



dN(t) 
dt 



-fN m (t) 



1 - 



N(t - t) 
K 



(10) 



which can also be treated as a generalization of the Hutchinson delayed logistic equation (jTl). 
This time-delay equation is capable of simulating a population that can essentially surmount 
the carrying capacity K for a limited period of time, then dropping below it subsequently [T9] . 
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(xi) Singular Malthus equations 

All equations enumerated above produce bounded solutions. In some cases, the pop- 
ulation dynamics seems to follow a law ending with divergent solutions. The well known 
Malthus equation [20] gives the exponential population growth. But sometimes, the popula- 
tion dynamics develops a super-exponential behavior, diverging at a finite time t c according 
to a power law of the time to the singular time t c . The simplest way of modeling such a 
behavior is by the equation 

dN 

with the power m > 1. This is a direct generalization of the Malthus equation, capturing the 
positive feedback of the population on the growth rate: the larger the population, the higher 
the growth rate. For m — 1, one recovers the usual Malthus equation with the exponential 
solution. When m > 1, the solution to the equation (TTTT) is of power law 



N(t) 



C 



(tc-t) Ve 



where 



The solution diverges hyperbolically at the critical time t c . Such strongly singular solutions 
were first discussed by von Foerster et al. [21] and applied to rationalize the super-exponential 
growth of the human world population [22TT25] . population dynamics and financial markets 
material failures and earthquakes [29, 30J, climate and environmental changes [3"TT 
and dynamics of other systems [5BT1S7] . In ecology, the correlation between population 
density and the per capita growth rate is known as the Allee effect [38J. The feedback 
between the population density, associated with the Allee effect, can lead to the increase of 
the effective growth rate, in the case of sufficiently large populations, or to the rate decrease 
and species extinction for small-density populations [5BTHTj . 

In addition to the differential equations describing population dynamics, there exist as 
well difference equations [121B3] and integro-differential equations [1HH5]. There are more 
complicated equations characterizing several factors, such as the population density, mass 
or weight dependence of individual members of different species, the dynamics of available 
food, and so on [IIJ[I5J[I7]. Some study the influence on population dynamics of available 
information [48J . It is also possible to investigate the spatial dependence of populations [1H] • 
More details on these and other types of equations can be found in the review articles [50] - f53] . 

We do not consider here the complications caused by the desire to take into account many 
various features of the studied populations. Our aim here is different: we concentrate our 
attention on the influence of the functional dependence of carrying capacities on population 
densities. Since this idea is rather new, it is necessary, first of all, to study the related 
effects for simpler equations, without overloading them by secondary specifications. Once 
the main influence of the functional carrying capacities is understood, it will be possible to 
complicate the equations by taking into account more and more mechanisms and specificities. 
For the same reason, the parameters characterizing the interactive species are treated as 
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fixed. In reality, the characteristics of each given biological species vary in general with the 
age of the individuals composing the group. However, it is always admissible to divide the 
populations into different age ranges characterized by similar birth and death rates. Another 
possibility is to consider each species being characterized by effective averaged parameters, 
which corresponds to what is called the mean-field approximation [36 l [37 j. l54] . 

In the dynamics of any population, one can distinguish several time scales. The shortest 
time scale is the interaction time tint, describing interactions between individuals. This time 
is much shorter than the observation time t jj S , during which the population is investigated. 
In order that the studied species could be characterized by fixed parameters, such as birth- 
death rates, the observation time should be shorter than the variation time t var during which 
individuals experience noticeable changes related, e.g., to their age. In this way, we keep in 
mind the situation, when 

Lint < tabs < 

Then all population parameters, except the carrying capacity, can be kept fixed. Our main 
point is that the carrying capacity can vary due to mutual interactions of individuals, hence 
it varies during the interaction time and this variation needs to be taken into account. 

The standard situation in treating the population evolution by means of the equations of 
the above types is that the carrying capacities are kept as fixed constant parameters. In the 
following Section 2, we propose an approach where the carrying capacities are functionals 
of the species populations. This makes it possible to drastically extend the applicability of 
the population-evolution equations to various situations describing the regimes that where 
unavailable with other models. The basic idea of the approach has been formulated in our 
previous papers [391IIQ], where some particular models were considered. Now, in Section 2, 
we propose a general framework for generating a large class of such models. In Sections 3 
and 4, we consider particular variants of the suggested equations, corresponding to those of 
Refs. [39JH0]. The difference from the previous works is three-fold. First, we simplify the 
consideration by a convenient choice of the scaling for the terms describing the carrying ca- 
pacities, which allows us to make a more straightforward classification of different dynamical 
regimes. Second, we emphasize the conditions under which solutions arise that are charac- 
terized by extreme events, such as the finite-time death and finite-time singularity of the 
species. Third, we suggest a new interpretation for an extreme event such as the finite-time 
singularity. The novel interpretation is based on the leverage effect and considers the sin- 
gularity as a manifestation of an evolutional boom followed by a crash. Our considerations 
are phrased for applications to the development, growth and possible collapse of biological 
as well as human societies, as they both follow similar dynamics with analogous underlying 
mechanisms [37] . 

2 Population evolution with functional carrying capac- 
ity 

The idea that the carrying capacity may be not a constant but a function of population 
fractions has repeatedly appeared in the literature in the form of general discussions. In Sec. 
2.1, we give a historical overview of these ideas that provide a firm justification for their 
mathematical representation in Sec. 2.2 and in the following sections. 
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2.1 General meaning of functional carrying capacity 

The carrying capacity of a biological species in an environment is generally understood as the 
maximum population size of the species that the environment can sustain indefinitely, given 
the food, habitat, water and other necessities available in the environment. In population 
biology, carrying capacity is defined as the environment maximal load [S3] , which is different 
from the concept of population equilibrium. Historically, carrying capacity has been treated 
as a given fixed value [56j|57]. But then, it has been understood that the carrying capacity 
of an environment may vary for different amounts of species and may change over time due 
to a variety of factors, including food availability, water supply, environmental conditions, 
living space, and population activity. 

Thus, the carrying capacity of a human society is influenced by the intensity of human 
activity, which depends on the level of technological development. When prehistoric hu- 
mans first discovered that crude tools and weapons allowed greater effectiveness in gathering 
wild foods and hunting animals, they effectively increased the carrying capacity of the en- 
vironment for their species. The subsequent development and improvement of agricultural 
systems has had a similar effect, as have discoveries in medicine and industrial technology. 
Clearly, the cultural and technological evolution of human socio-technological systems has 
allowed enormous increases to be achieved in carrying capacity for our species. This in- 
creased effectiveness of environmental exploitation has allowed a tremendous multiplication 
of the human population to occur [58||59]. 

Technology is an important factor in the dynamics of carrying capacity. For example, 
the Neolithic revolution increased the carrying capacity of the world relative to humans 
through the invention of agriculture. Currently, the use of fossil fuels has artificially increased 
the carrying capacity of the world by the use of stored sunlight, albeit with increasingly 
negative externalities, such as global warming, ocean acidification and the indirect reduction 
of diversity. Other technological advances that have increased the carrying capacity of the 
world relative to humans are: polders, fertilizer, composting, greenhouses, land reclamation, 
and fish farming. Agricultural capability on Earth expanded in the last quarter of the 20th 
century. Whether this is sustainable is debatable. There are signs that human-induced soil 
erosion as well as destabilization of sensitive ecosystems may lead, at the same time, to a 
reduction of agricultural capability over the coming decades, such as in Africa where the 
population is expected to double before 2050. The change in the carrying capacity of the 
habitat and environment supporting a human society can be described by the consumption 
impact [60J, which is proportional to the population size, with a coefficient characterizing 
the technology level. 

One way to estimate human influence on the carrying capacity of the ecosystem is to use 
the so-called ecological footprint accounting method that provides empirical, non-speculative 
assessments of human activities with regard to the preservation or destruction of the Earth 
carrying capacity. It compares regeneration rates (biocapacity) against human demand (eco- 
logical footprint) in the same year. The results show that, in recent years, humanity demand 
has exceeded the planet biocapacity by more than 20 percent [51]. The present situation of 
rapid population growth in some regions, massive overexploitation of resources and steady 
accumulation of pollution and wastes diminishes the Earth carrying capacity. To a first 
approximation, with all the caveats associated with the heterogeneity of technological devel- 
opments in different parts of the World, one can consider an average footprint per person, 
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which leads to an estimation of the decrease of the Earth carrying capacity as roughly pro- 
portional to the size of the total population. The question is how and by what means this 
change of the Earth carrying capacity for humanity will pay back [62] . 

Mutual coexistence and symbiosis of several species also strongly influences the carry- 
ing capacities of the species, with the changes being, to a first approximation, proportional 
to the species numbers. For example, humans have increased the carrying capacity of the 
environment for a few other species, including those with which we live in a mutually ben- 
eficial symbiosis. Those companion species include more than about 20 billion domestic 
animals such as cows, horses, pigs, sheep, goats, dogs, cats, and chickens, as well as certain 
plants such as wheat, rice, barley, maize, tomato, and cabbage. Clearly, humans and their 
selected companions have benefited greatly through active management of mutual carrying 
capacities [63] . 

Interactions between two or more biological species are known to essentially influence 
the carrying capacity of each other, by either increasing it, when species derive a mutual 
benefit, or decreasing it, when their interactions are antagonistic [6H466] . The same applies to 
economic and financial interactions between firms, which also form a kind of symbiosis, where 
the interacting firms develop the carrying capacity of each other also roughly proportionally 
to their sizes [67] . 

Many authors (e.g., Del Monte-Luna et al. [68]) stress that, due to the influence on the 
carrying capacity resulting from the existing populations, its original definition implying a 
constant value has lost its meaning. As a consequence of the feedback loops of the population 
sizes, the notion of carrying capacities has taken a broader sense. Carrying capacity should 
be understood as a non-equilibrium relationship or function that depends on the population 
size and on the symbiotic relations between the interacting populations. It characterizes the 
growth or development of available resources at all hierarchical levels of biological integration, 
beginning with the populations, and shaped by processes and interdependent relationships 
between finite resources and the consumers of those resources [65] . 

The above discussion shows that, in general, the carrying capacity is not a fixed quantity, 
but it should be considered a function of population sizes. In the case of a single species, 
the carrying capacity can be created or destroyed by the species activity. The simplest as- 
sumption is to take the species impact as proportional to the species population size. When 
species coexist, their carrying capacities are influenced by the species mutual interactions, 
either facilitating the capacity development or damaging it. Being functions of the species 
populations, such nonequilibrium carrying capacities can be naturally represented as poly- 
nomials over the population numbers [69J. 

2.2 Mathematical formulation of basic equations 

Let the considered society consist of several species enumerated by an index i = 1,2, . . .. 
The number of members in a society is denoted as iVj = Ni(t), 

The main idea of our approach is that the carrying capacity of each species is not a fixed 
constant, but it is a functional 

Ki = Ki({Ni}) (12) 

of the set {Ni} of the species populations. This assumption takes into account that the 
species may interact not merely directly but also by influencing the carrying capacities of 
each other as well as their own carrying capacity. 
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A general form of the evolution equations, that takes into account both direct interactions 
and mutual influence on the carrying capacities, can be written as 



dt 



N t , (13) 



where is the functional carrying capacity f)12p . If the latter were a fixed parameter, 
one would return to the generalized Lotka-Volterra predator-prey model ([3]). But in our 
approach, the carrying capacities are nontrivial functionals of populations. 

In a particular case, when the species do not display strong direct interactions, but their 
mutual correlations are mainly through influencing the carrying capacities of each other, 
then Eq. ([TBI reduces to the evolution equation 

dt h Ki y ! 

Here the effective rate 

is the difference between the birth and death rates of the corresponding species. When birth 
prevails over death, then 7, > 0, while if death is prevailing over birth, then 73 < 0. In 
economic applications, birth translates into gain and death into loss. The parameter 

d = cf mp - cr p (i6) 

characterizes the difference between the competition and cooperation of the members inside 
a given type of species. When competition is stronger than cooperation, then C{ > 0, while 
if cooperation is stronger, then d < 0. 

It is important to stress that Eq. ({TBI) is principally different, both mathematically as 
well as by its meaning, from the predator-prey model Similarly, Eq. (j!4p is principally 
different from the logistic equation §6§ or its variants ([7]) and (JHD- As a result, the equations 
with functional carrying capacities can display novel types of solutions allowing for the 
consideration of effects that are absent in other evolution equations. In the following sections, 
we consider concrete examples of evolution equations with functional carrying capacities. 



3 Action of society on its own carrying capacity 
3.1 Justification for equation form 

In the previously studied variants of the logistic equation, the carrying capacity is treated 
as a given quantity. However, it is often the case that the society activity does influence 
its own carrying capacity that can be either enhanced by producing new goods, materials, 
knowledge, and so on, or can be destroyed by unreasonable exploitation of resources, e.g., by 
deforestation, polluting water, and spoiling climate. Therefore, the carrying capacity, taking 
into account such feedback effects, must be a functional K = K(N) of the population N. 
Thus, the population evolution is characterized by the equation 

dN _ CN 2 

~dt~ 1 ~ kJn) ' 
9 



(17) 



with the carrying capacity depending on N itself. 

Moreover, the creation or destruction of the carrying capacity by the society members 
does not occur immediately, but is delayed since any creation or destruction requires time 
for its realization. Hence, the variable N, entering the carrying capacity, should be delayed 
in time by a lag r, so that K(N) = K(N(t — r)). 

Different types of carrying capacity can be introduced that depend on a delayed popula- 
tion variable. In the present paper, we consider a simple linear form 

K(N) = A + BN(t - r) . (18) 

Here the first term A > is a natural carrying capacity, provided by Nature. The second 
term is the created or destroyed capacity, depending on whether the society activity is 
constructive or destructive. The parameter B is the production factor, if it is positive, and 
it is a destruction factor, when it is negative. Form ffTB"]) of the carrying capacity agrees 
with the assumption of additivity, when its different parts sum to produce the total carrying 
capacity. 



3.2 Reduced quantities and choice of scaling 

As the population numbers can be very large, it is therefore more convenient to deal with 
the reduced quantities 

N(t) , , 

(19) 



x 



X, 



eff 



measured in units of some typical population size N e ff. It is also convenient to introduce 
dimensionless parameters for the natural carrying capacity 



.4 



N, 



eff 



c 



and for the production-destruction factor 

b = B 

The total dimensionless carrying capacity 



C 



y 



K(N) 



N, 



eff 



c 



(20) 

(21) 

(22) 
(23) 



takes the form 

y = a + bx(t — t) . 

Up to now, the effective value N e ff has been arbitrary. By a special choice of the scaling, 
it is possible to simplify the equations and to make a more transparent classification of 
arising dynamical regimes. It is convenient to choose 



Neff = A 

Then parameters (I2"U1) and (I2T1) reduce to 

a = 1 , b - 



C 



B 
A 



N, 



eff 



(24) 



(25) 
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The dimensionless carrying capacity 023 j) becomes 

y = 1 + 6z(t - r) . (26) 

Let us also define the signs of the birth-death rate and of the competition-cooperation 
parameter as 

7 C 
o-i = sgn7 = r-r , <r 2 = sgnC = — . (27) 

n\ \ c \ 

Depending on these signs, the following situations can occur. 

(gain + competition) , 

(gain + cooperation) , , . 



o x = 


+1, 


2 = 


+ 1 


0\ = 


+1 , 


02 = 


-1 


a i = 


-1 , 


2 = 


+ 1 


Oi = 


-1, 


2 = 


-1 



Using the above notations and measuring time t > in units of I/7, we reduce Eq. (|T7|) 
to 2 

— = 01X - 2 — , (29) 
at y 

with the carrying capacity ( |26|) . This equation is to be complemented by the initial conditions 

x(t) =x (t < 0) , 

y(*) = 1/0 = 1 + bx (t<0). (30) 

The solution for x, by its meaning, is to be positive. The production-destruction factor b 
can take any real values, being positive for the constructive society activity, while negative, 
for its destructive activity. 

3.3 Evolutionary stable states 

One of the most important problems in studying any evolutional model is the determination 
of evolutionary stable states. These are given by the stable stationary solutions to the 
considered equation. In order to analyze the stability of the solutions to the differential delay 
equations, we employ the Lyapunov stability theory following the work by Pontryagin [70] 
and the books mmaHSI. 

The stationary states of Eq. (|2~91 are defined as the solutions to the fixed-point equation 

2 (x*) 2 . , 

o-ix* = . 31 

1 + ox* 

This yields two fixed points 

x\ = , x\ = ^— . (32) 
02 - b 

Resorting to the Lyapunov stability analysis, we need to consider a small deviation 

8xj(t) = xj(t) - x* (j = 1, 2) (33) 
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from the related fixed point. This deviation satisfies the equation 



j 6xj(t) = CjSxjit) + DjSxjit - r) , (34) 



in which 



2a 2 x* 



X; 



a =<Ti - — , Dj = ba 2 1 



1 + bx* J 1 + bx* 



For the corresponding fixed points, these parameters are 

C 1 = a 1 , D x = , 
6(oi - 1) - cr 2 6(T 2 

^2=CTlT7 TT - > -^2 



6(0-! - 1) + a 2 ' z [6(0-! - 1) + a 2 ] 2 ■ 

Looking for the deviation in the exponential form 

Sxj(t) oc e Xjt , 

we obtain the equation 

A, = Cj + D je -^ (35) 
for the characteristic exponents Xj. By using the notation 

Wj = (A, - Cj)t , Zi ee r^e-^ T , 

equation (1351) becomes 

which is nothing but the equation defining the Lambert function Wj. Therefore, the characteristic- 
exponent equation (|35|) acquires the form 

\ J = C j + ±W j (D J e-°> r ). (36) 



The stationary solution is stable when the real part of the characteristic exponent is negative, 
< 0. 

To proceed further, we shall analyze separately the cases listed in Eq. ( 1281) . 



4 Society with gain and competition 

Under the prevailing gain (birth) and competition, when 

cri = l, o- 2 = l, (37) 
the evolution equation ( 1291) reads as 

- x(t ) _ ^ ( 38 ) 

dt ~ X{t) 1 + bx(t - r) ■ [6H) 



12 



At the initial stage, when < t < r, we have the exact solution 

so(l + te )e* 
* W l + Xo ( & -l + e *) l"St<r;. 

This can be used for constructing by iteration an approximate solution at the second stage, 
when r < t < 2t. Then we could find an approximate solution at the third step, and so on. 
However, the accuracy of such iterative constructions quickly deteriorates and is admissible 
only over a couple of initial steps. More accurate solutions are to be found by numerically 
solving the evolution equation (!38il . 

The stability analysis of the previous section shows that the fixed point x\ = is unstable 
for all b and r. The second fixed point 



is stable when either 
or when 
where 

Tn = 

The stability region is shown in Fig. 1. 

Varying the system parameters and initial conditions, we can meet the following dynamic 
regimes. 



- 1 < b < 1 , r > , (40) 
b < -1 , t < T , (41) 



t = —^=^= arccos ( - ) . (42) 



4.1 Punctuated unlimited growth 

For the parameters 

b > 1 , t > (x > 0) , (43) 

the population grows by steps, as shown in Fig. 2. The growth continues to infinite times. 
This is a typical example of the punctuated evolution caused by the fact that the production 
factor b is positive and sufficiently large. Hence, the carrying capacity is produced by the 
population, with a delay r. 



4.2 Punctuated growth to stationary state 

The punctuated growth is not always unbounded, but it can be bounded by the fixed point, 
provided the initial condition xq is smaller than x* and the parameters are 

< b < 1 , r > (x < x*) . (44) 

This regime is presented in Fig. 3. 
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4.3 Punctuated decay to stationary state 



When the parameters are the same as in Eq. fl44j) . but the initial condition xq is larger than 
the fixed point x*, 

x > x* = - , (45) 

l — o 

then there appears the punctuated decay, as illustrated in Fig. 4. 

4.4 Punctuated alternation to stationary state 

If the carrying capacity is destroyed by the population, then there can occur a punctuated 
alternation to a stationary state, when the parameters and the initial condition are 

-1<6<0, r>0 r 0< |^f) ■ (46) 
This is depicted in Fig. 5. 

4.5 Oscillatory approach to stationary state 

For sufficiently large destruction factor, there arises a regime of an oscillatory approach to 
a stationary state, as is presented in Fig. 6. This happens under the parameters and the 
initial condition being defined by the inequalities 

b < -1 , r < r (x < y^iJ , (47) 



where the lag r is defined in Eq. (|42|) . 

4.6 Sustained oscillations 

There exists a lag t\ = ti(o), such that, when 

6<-l, ro<T<n ( Xo< |^) ' ( 48 ) 

then oscillations do not decay, but continue without attenuation, as in Fig. 7. The lag T\(b) 
can be found only numerically. 

4.7 Punctuated alternation to finite-time death 

If the time lag surpasses the value T\ = Tx(6), the alternating solution exists only for a limited 
time. At the death time ta, given by the equation 

1 + bx(t d - r) = , (49) 

all population becomes extinct, as in Fig. 8. This happens for the parameters 

b < -1 , t > n (x < y^iJ . (50) 
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4.8 Growth to finite-time singularity 

In the case, where the activity of the population is destructive, time lags are large, and the 
initial condition is also large, so that 

b < , t>t c (x > ^| J , (51) 

the population dynamics becomes dramatic, diverging at a finite time, called the critical 
time. The divergence is hyperbolic, according to the law 

x(t) oc — — - (t -> t c - 0) , (52) 
t c t 

as is demonstrated in Fig. 9. The values of the critical lag r c and the critical divergence 
time t c can be found numerically. 



4.9 Unlimited exponential growth 

For shorter time lags, when 

b < , r<r c L 0> ±-\ (53) 

the divergence moves to infinity, the solution being a simple growing exponential, as shown 
in Fig. 10. 

In this system with gain and competition, there may happen two extreme events, the 
finite-time death at a death time td and the finite-time singularity at a critical time t c . These 
two extreme events occur under the condition of a destructive activity of the population. 
The finite-time death is caused by the destruction of all resources. The finite time singularity 
implies that close to this critical point, the dynamic regime has to be changed, according to 
the accepted interpretation of such singularities [27JI321I1D]- Such a change of the dynamic 
regime is analogous to the occurrence of critical phenomena in statistical systems [36j[37J 
EUEEj. An interpretation of the finite-time singularity, based on the leverage effect, will be 
given below. 



5 Society with gain and cooperation 

When gain prevails over loss, and cooperation over competition, that is, when 

ax = 1 , (7 2 = -1 , (54) 

the evolution equation takes the form 

dx(t) . . x 2 (t) 

= x (t) + r • 55 

dt w 1 + bx(t - t) y ' 

There are no stable stationary solutions in that case. Depending on the system parame- 
ters and initial conditions, there can arise the following dynamic regimes. 
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5.1 Growth to finite-time singularity 

When the population activity is productive, but the time lag is long, so that 

b > , t > t c (x Q > 0) , (56) 
or if the activity is destructive, when 

b < , r > (x < ttj ) , (57) 



|6| 

then the solution diverges at a finite critical time t c . The behavior is the same as in Fig. 9. 

5.2 Unlimited exponential growth 

Productive activity, under cooperation and not too long time lags, such that 

b > < t < t c (x Q >0), (58) 
result in an exponential growth, as in Fig. 10. 

5.3 Punctuated unlimited growth 

For the parameters 

6<-l, r>0 ^o> j^prr) > ( 59 ) 
the solution displays unlimited punctuated growth, as in Fig. 2. 

5.4 Punctuated decay to finite-time death 

Destructive activity, under one of the conditions, when either 



- 1< b < , r > (^o > Y^\b\) ' (60) 

or when 

6<-l, r>0 ( x o<j^ =r j^ , (61) 

leads to population extinction, at the death time given by Eq. (H9|) . But the dynamics 
for this CclSG, clS IS shown in Fig 11, is different from that of Fig. 8. In the present case, 
there are no alternations, but the decay to zero is monotonic, exhibiting a finite number of 
quasi-plateaus. 

Under the conditions of gain and cooperation, there are two types of extreme events, 
the finite-time singularity at a critical time t c and the finite-time death at a death time td- 
The finite-time death is caused by the destructive population activity. And the finite-time 
singularity means that, close to the singularity point, the system experiences a change of 
dynamic regime. 
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6 Society with loss and competition 

Under prevailing loss and competition, when 

a x = -1 , a 2 = 1 , (62) 

the evolution equation becomes 

^& - _ x(t) _ x " (*) f63 ) 
dt ~ 1 J 1 + 6x(t - r) • [b6) 

There are two stable fixed points. One is the trivial point 

x{ = 0, (64) 

which is stable for all parameters 

- oo < 6 < oo , r > . (65) 

Another stationary point 

T* = - 

1 + b 



x* 2 = --^— = x* (66) 



is stable for the parameters 

6<-l, r<r , (67) 

where 

T ° s "«» (]if) • < 68 > 

Thus, there is the bistability region shown in Fig. 12. Solutions tend to one of the two sta- 
tionary states, when the initial conditions are in the basin of attraction of the corresponding 
fixed point. The following regimes can arise. 

6.1 Monotonic decay to zero 

In a society with prevailing loss and competition, the decay to zero, as in Fig. 13, seems to 
be a natural type of behavior. This happens when either 

b > , t > (x > 0) , (69) 

or when 

b > , r > U<ij . (70) 

6.2 Oscillatory convergence to stationary state 

When the parameters are such that 

b < -1 , < r < r (x > jj^—^J , ( 71 ) 

the population fraction x oscillates in time, converging to the stationary state ( |66|) . as is 
shown in Fig. 14. Oscillations are caused by the presence of the time delay. 
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6.3 Everlasting nondecaying oscillations 

For the parameters 

b < -1 , r < r < n hr > j^py) > ( 72 ) 

the solution oscillates without decay, similarly to the behavior in Fig. 7. The time lag tq is 
given by Eq. (1681 and ri is defined numerically 

6.4 Punctuated growth to finite-time singularity 

A rather interesting behavior of the population dynamics happens for the parameters 

b < -1 , n<T<T 2 (x > Jj^—^) ■ (73) 

Then the solution experiences several punctuations, after which it diverges, as is illustrated 
in Fig. 15, at the critical time defined by the equation 

1 + bx(t c - r) = . (74) 

When the final rise is preceded by a fall, this behavior is reminiscent of the Parrondo effect 



6.5 Up-down convergence to stationary state 

A highly non-monotonic behavior exists for the parameters 

6<-l, 0<t<t c (-^ < x < j^z^j > (75) 

where the time lag r c can be found only numerically. In this case, the solution, first, bursts 
out upwards, after which it decays to the stationary value x*, as in Fig. 16. 

6.6 Growth to finite-time singularity 

Under the parameters 

K-l, r>r c (L <X0<] J-y (76) 

the solution diverges at a finite critical time, without any punctuation, in the same way as 
in Fig. 9. 
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6.7 Unlimited exponential growth 

In the region of the parameters 

-1<6<0, 0<r<r c (x > , (77) 

the solution grows exponentially, as in Fig. 10. 

For a society with prevailing loss and competition, there are two extreme events, both 
characterized by a finite-time singularity at a critical time t c . These regimes occur under a 
strong destructive activity of the population and a rather long time lag. The divergence can 
be understood as a critical point where the society dynamics qualitatively changes. 



7 Society with loss and cooperation 

When loss and cooperation prevail, so that 

<7=-l, (T 2 = -I, (78) 

the population evolution equation is 

d J^l = _,(() + fW . (79) 
dt w 1 + bx(t - t) y 1 

There exists the sole evolutionary stable state 

x* = (80) 

that is stable for all parameters 

-oo<6<oo, r>0. (81) 
The following dynamic regimes are possible. 

7.1 Monotonic decay to zero 

For the initial conditions in the attraction basin of the stable fixed point, the solutions decay 
to zero with time, as in Fig. 13. This happens when either 

b<0, r>0 ^ 0<r -L_) , (82) 

or when 

b > 1 , t > (x > . (83) 
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7.2 Growth to finite-time singularity 



If the initial conditions are outside of the attraction basin of the fixed point flHUj) . they can 
diverge at a finite critical time, similarly to the behavior in Fig. 9. This happens when 
either 

b<0, r>0 ( Xo< ^) ' (84) 

or when 

< b < 1 , r>r c (x Q > ——r ] . (85) 



1-b 



7.3 Unlimited exponential growth 

For the parameters 



C)</>< I . r t, [ Xo> -L-] . (86) 



the solution exhibits exponential growth, as in Fig. 10. 



7.4 Monotonic decay to finite-time death 

Finally, for the parameters 

6<0, 0<r<r d (x 0> ±-J , (87) 

the population becomes extinct at a finite death time, as in Fig. 17. The death time is 
defined by an equation having the same form as Eq. (H9|) . However the decay to death now 
is monotonic, which distinguishes it from the punctuated behavior before death, shown in 
Fig. 8 and Fig. 11. 

The society with prevailing loss and cooperation can exhibit the finite-time singularity as 
well as the finite-time death. These two types of extreme events happen under the destructive 
activity of population. 

Summarizing, all extreme events, except one, occur when the population destroys its 
carrying capacity. The sole exception is the case of a society with gain and cooperation, when 
there can arise a finite-time singularity under b > 0, i.e., when the activity of the population 
is productive. This latter type of finite-time singularity is analogous to that studied in 
Ref. |27] . Its appearance means that, near the critical time, the society becomes unstable 
and requires to change its parameters, for instance replacing cooperation by competition. 
It seems to be rather clear that, when the population grows too much, the competition of 
individuals must come into play, becoming prevailing over their cooperation. The finite-time 
singularities, occurring under the destructive society activity, imply the existence of some 
critical events, whose detailed interpretation will be given in Sec. 11. 
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8 Mutual influence of symbiotic species on their car- 
rying capacities 



8.1 Classification of symbiosis types 

When the considered society is structured with several species, it is necessary to character- 
ize their interactions. The standard way of doing this is by assuming the equations of the 
predator-prey type (j3J), with direct interactions of species that eat each other. Such equa- 
tions, however, cannot describe indirect interactions, when the species do not kill each other, 
but influence the carrying capacities of each other. Therefore, the predator-prey equations 
are suitable for describing the predator-prey relations, but are not suitable for characterizing 
symbiotic relations [3D] . 

Examples of symbiosis are ubiquitous in biology and ecology [731476] . It is also widespread 
in human societies. For example, one can treat as symbiotic the interrelations between firms 
and banks, between population and government, between culture and language, between 
economics and arts, and between basic science and applied research. 

Considering purely symbiotic relations, we need equation ( I14p . with the carrying capaci- 
ties being functionals of the species populations. The natural form of such carrying capacities 
for symbiotic species is 

Ki = Ai + BiSidNi}) . (88) 

Here Ai > is the natural carrying capacity, provided by nature, for an z-th species. The 
coefficient characterizes the strength of influence of other species on the carrying capacity 
of the z-th species. When Bi is positive, it can be called the production factor, while, if it is 
negative, it is the destruction factor. The function S({Ni}) is a symbiotic function specifying 
the mutual relations between symbiotic species. Since the sign has already been attributed 
to the factor Bi, the symbiotic function can be treated as non-negative. 

Depending on the kinds of symbiotic relations, that is, on the signs of the factors Bi, 
there can occur different variants of symbiosis. To illustrate this, let us analyze the case of 
two symbiotic species for which there can exist the following types of symbiosis. 

(i) Mutualism, when both species are useful for each other, developing their mutual 
carrying capacities: 

B 1 > , B 2 > (mutualism) . (89) 

(ii) Parasitism, when one of the species is harmful for another, or both species are 
harmful for each other, destroying the carrying capacities, which happens under one of the 
pairs of inequalities below: 

B x > , B 2 < , 

B x < , B 2 > , (parasitism) (90) 

B x < , B 2 < . 

(iii) Commensalism, when one of the species is useful for another, while the latter is 
indifferent to the existence of the first species, which corresponds to the validity of one of 
the pairs of equations: 

Bx > , B 2 = , 

B x = , B 2 > (commensalism) . 
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8.2 Normalized species fractions 

We continue analyzing the symbiotic coexistence of two kinds of species. As always, it is 
more convenient to work with reduced quantities. So, we introduce the reduced fractions 

% = Tr , Z = y— , (92) 

iV e// ^eff 

whose normalization values N e ff and Z e ff will be chosen later. We define the dimensionless 
carrying capacities 

Vl = n at ' V2 = 7^7 — ( 93 ) 

and the relative birth rate 

a = — . (94) 
7i 



With these notations, the symbiotic equations f)T4j) . in the case of two types of species, 
reduce to 

dx x 2 dz z 2 

— = x , — = az , (95) 

dt yi dt y 2 

where time is measured in units of I/71. By their definition, the solutions x and y are 
non-negative. The equations are complemented by the initial conditions 

ar(0) = x , z{0) = z . (96) 

For the following analysis, it is necessary to make concrete the explicit forms of the 
carrying capacities flHHj) . 



9 Symbiosis with mutual interactions 
9.1 Derivation of normalized equations 

The action of the species on the carrying capacities of each other can be different, depending 
on whether, influencing the carrying capacities, the species interact or not. If the species, 
in the process of influencing their carrying capacities, interact with each other, then the 
carrying capacities (188]) can be represented in the form 

K X =A X + B X N X N 2 , K 2 = A 2 + B 2 N 2 N X . (97) 

Generally, the populations N x and A^ 2 in these carrying capacities could depend on the shifted 
time, when one would have 

Kt = Ai + BiN^t - Ti)Nj(t - T j) , 

where i ^ j. However, we need, first, to understand the influence of symbiosis without the 
time lag. Therefore, we consider below the interactions without time delay. 
Introducing the dimensionless natural carrying capacities 
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and dimensionless symbiotic factors 

^Ti^/, g^Jl^kl (99) 

translates Eqs. ( 197]) into the dimensionless expressions 

2/i = fl i + bxz , y 2 = a 2 + gxz . (100) 

Since the scaling values N e ff and Z e ff are arbitrary, it is reasonable to choose them so 
as to simplify the equations. For this purpose, we set 

AW, - ^ , - ^ ■ (101) 

Then the natural carrying capacities ( 198|) become 

a 1 = a 2 = l . (102) 
And the total carrying capacities f llOOp read as 

yi — 1 + bxz , 1/2 = 1 + gxz . (103) 

As usual, we measure time in units of I/71. The most interesting case in symbiosis is 
when the species influence each other throughout their lifetimes, and when these lifetimes 
are of comparable durations. If this were not the case, i.e., with very different lifetimes, the 
symbiotic relations could not be supported for a duration longer than the shortest lifespan, 
making symbiosis inefficient for the longer-lived species. Therefore, we assume that the 
symbiotic species have comparable growth rates, because the inverse of the growth rate sets 
the time scale of lifetime, and the later is often found proportional to the growth period, at 
least for mammals [77]. We thus set a = 1 and obtain the equations 

dx x 2 dz z 2 . 

-— = x - - — , — = z - - . 104 

at 1 + bxz dt 1 + gxz 

We can note that these equations are symmetric with respect to the simultaneous interchange 
between x with z and between b with g. This symmetry will result in the corresponding 
symmetry of the following solutions. 

9.2 Evolutionary stable states 

Again we use the Lyapunov stability analysis [IIJ[l6j[70] . Equations (1 104ft possess the non- 
zero stationary state 



* 1 
x = — 

2g 



l-b + g- ^(1 + b - g) 2 - 4b 



z * = ~k V + b ~ 9 ~ V^ + b-g) 2 -^] . (105) 

It is stable when either 

b < , -00 < g < +00 , (106) 
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or when 

< b < 1 , g < g c , (107) 

or when 

b > 1 , 5 < , (108) 

where the critical value g c is 

g c = (i-VP) 2 <i. (109) 

The stability region is depicted in Fig. 18. The basin of attraction of this stationary state, 
depending on the signs of the symbiotic factors, is defined by the following equations: 

x zo < (b < 0, g > 0) , 

XoZo < -rr {b > 0, g < 0) , 
\9\ 

x z < "wijj^i , jjjj (&<0, #<0). (110) 

The solutions to the symbiotic equations (I104p should be compared to those of the un- 
coupled equations 

^ = x-a; 2 , % = z~z 2 (b = 9 = 0), (111) 
ere etc 

corresponding to the case of no symbiosis, when the stationary states are x* — z* — 1. 

Solving numerically the system of equations (j!04p for different symbiotic factors and 
initial conditions yields the following possible dynamic regimes. 



9.3 Convergence to stationary states 

For the system parameters in the region of stability and for the initial conditions in the 
basin of attraction of the non-zero fixed point (11051) . both species develop and converge to 
the stationary state. This is illustrated in Fig. 19 for different types of symbiosis, where, for 
comparison, the solutions for the case of no symbiosis are also presented. The four possible 
cases are illustrated in the four panels of Fig. 19, depending on the relative positions of x(t) 
and z(t) compared with the solution of the uncoupled equations (II 111) . 

9.4 Unlimited exponential growth 

When stationary solutions do not exist, so that either 

< b < 1 , g > g c , (112) 

or when 

6>1, g>0, (113) 

or when they exist, but the initial conditions are taken outside of the attraction basin, then 
the populations of both species grow exponentially with time. 
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9.5 Finite-time death and singularity 



In the case of mutual parasitism, there can happen an extreme solution when one of the 
species becomes extinct at a finite critical time, while the other species displays a finite-time 
singularity, as is shown in Fig. 20. This happens when the initial conditions are outside of 
the attraction basin so that either 



\b\ 

or if 



1 < x z < — (b<g<0) , (114) 



t^t < x z < -j- (g < b < 0) . (115) 

\g\ V>\ 

The critical time is defined by one of the corresponding equations: 

x(t c )z(t c ) = -J- (b < g < 0) , 
\b\ 

x(t c )z(t c ) = — {g < b < 0) . (116) 

\g\ 

The appearance of such an extreme solution is caused by the mutual parasitism of the 
species, destroying the carrying capacities of each other. 



10 Symbiosis without direct interactions 
10.1 Derivation of symbiotic equations 

In many cases, symbiotic species influence each other by increasing (improving) the carrying 
capacities of each other, which does not involve direct interactions between the species. The 
most known example of this type is the symbiosis between tree roots and fungi. In that case, 
the carrying capacities (1551) can be written in the form 



K X = A X + B X N 2 , K 2 = A 2 + B2N, . (117) 

The dimensionless carrying capacities ( 1931) now read as 

y 1 = a 1 + bz, V2 = a 2 + gx. (118) 

Employing the scaling of Eqs. fl 1 1 1) gives normalization ( I102p and the carrying capacities 
( 11181) become 

Vl = l + bz, y 2 = l+gx. (119) 
Thus, we come to the symbiotic equations in dimensionless form 

dec dz z , x 

-r = x ~ --, — r- » -r = z ~ ; ■ 120 

dt 1 + bz dt 1 + gx 

There exists again the symmetry with respect to the simultaneous interchange between x 
and z and between b and g. 
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10.2 Evolutionary stable states 

Equations (11201) possess a non-zero stationary state 



which is stable when either 

or when 

where 



* 1 + b * 1 + 9 

1-05 1 - ^ 



- 1 < 6 < , g>-l, (122) 
& > , < g < g c , (123) 



g c = I ■ (124) 



The stability region is presented in Fig. 21. 

If the symbiotic relations correspond to mutualism or commensalism, then the attraction 
basin of the stationary solution (I12ip is the whole region of positive initial conditions: 

x > , z > (b > 0, < g < g c ) . (125) 

But if at least one of the species is parasitic, then the attraction basins are defined by one 
of the conditions, depending on the signs of the symbiotic factors: 

x < , z > (b > 0, g < 0) , 

\g\ 

x > , z < (6 < 0, g > 0) , 

X < 7-r , ^0 < 77T (6 < 0, 5 < 0) . (126) 
Ifl'l |o| 

The following dynamic regimes are possible. 

10.3 Convergence to stationary states 

If initial conditions are in the attraction basin, then both species converge to their stationary 
populations. The convergence can be monotonic or not, depending on the system parameters 
and initial conditions, as is demonstrated in Fig. 22. 

10.4 Unlimited exponential growth 

For the parameters outside the stability region, such that 

b > , g > g c , (127) 
there exists a solution with exponential growth in time for both species. 



26 



10.5 Finite-time divergence 

Extreme solutions appear when at least one of the species is parasitic and initial conditions 
are outside of the attraction basin. Thus, when either 

x > t^t , z > (6 > 0, g < 0) , (128) 

\g\ 

or when 

x > , z > i (6 < 0, (/ > 0) , (129) 

then one of the species experiences a finite-time singularity at a critical time t c that is defined 
numerically. In this case, when approaching t c , one of the following behaviors arise: 

x(t) —> x(t c ) < oo , z(t) — > oo , 

x{t) -> oo , z{t) -> z(t c ) < oo , (130) 

where the first line corresponds to conditions (|128p . while the second line, to conditions 
(11291) . The typical behavior of populations is shown in Fig. 23. 

10.6 Finite-time extinction 

Parasitic symbiotic relations may end with one of the species being extinct and the other 
continuing its life without symbiosis. When either 

x < , z > (b > 0, g < -1) , (131) 

\g\ 

or when 

x < , > ]i (6 < 0, 5 < 0) , (132) 

\g\ \ b \ 

then the species z dies at a finite time tj, defined by the relation 

x{t d ) = ir . (133) 

That is, the species x kills the species z: 

x(t) -> s^) , z(t) ^0 (t -)• t d ) . (134) 
The opposite situation, when the species z kills the species x occurs if either 

x > , z < (6 < -1, p > 0) , (135) 

or if 

x > T^T , < ^ (6 < 0, g < 0) . (136) 

\g\ \b\ 

Then the species x dies at a finite time given by the relation 

<U) = ||| , (137) 

so that 

x(t) -> , z(t) -)• z(t d ) (t -> t d ) . (138) 
The corresponding behavior is illustrated in Fig. 24. 
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11 Interpretation of extreme events in population evo- 
lution 



11.1 Types of extreme events 

In the population evolution, there may happen two types of extreme events, finite-time death 
and finite-time singularity. The origin for the occurrence of finite-time death is rather clear. 
This happens when the carrying capacity of the species is destroyed, either by the species 
themselves or by the parasitic symbiosis of other species. The destroyed carrying capacity 
makes it impossible the long-term existence of the species that, thus, go towards extinction. 

Finite-time singularity can be due to two causes. One reason is the existence of cooper- 
ation between the members of species, as in Sees. 5.1 and 7.2. This type of the finite-time 
singularity means that the society, in which cooperation persists under fast increasing num- 
bers of its members, becomes unstable and, to be stabilized, requires that cooperation be 
changed into competition. The necessity for such a change looks rather evident and is easily 
understandable. Really, in the presence of a strongly increasing population, the competition 
for the means of survival will become unavoidable. 

A more elaborate mechanism operates in the case of the finite-time singularities occurring 
under competition, as in Sees. 4.8, 6.4, 6.6, 9.5, and 10.5. In all these cases, the singularities 
appear under the destruction of the carrying capacities either by the society itself or by 
a parasitic symbiotic species. It may seem quite strange that, while the carrying capacity 
is being destroyed, the population continues growing. To understand the origin of such a 
paradoxical effect and of these finite-time singularities, let us consider in turn the different 
types of finite-time singularities found in Sees. 4.8, 6.4, 6.6, 9.5, and 10.5. 

11.2 Boom and crash in society with gain and competition 

This corresponds to the case studied in Sec. 4.8, of a finite-time singularity occurring at a 
critical time t c . The divergence is of the hyperbolic type (|52|) . This extreme event happens 
under the destructive action of the society on its own carrying capacity, when b < 0. The 
parameters are such that, at the initial moment of time, the effective carrying capacity is 
negative, 

2/(0) = 1 - \b\x < . (139) 

How would it be possible to understand the existence of a negative carrying capacity? 
For some simple biological species, as ants or bees, the negative capacity would, probably, be 
impossible. Such species would not be able to live at all. However, for more complex societies 
such as human societies, the negative carrying capacity may have sense. For instance, 
humans do extract non-renewable resources that become progressively exhausted forever, 
they destroy their habitat, poison rivers, pollute air, cut forests, and so on. At the same 
time, humans possess the ability of regenerating the habitat by cleaning rivers, or even 
oceans, and planting trees. Thus, humans may spoil their habitat to such an extent that it 
would require a hard work for its recovering. In that sense the effective carrying capacity 
can become negative for a while, implying the necessity of its recuperation in the positive 
domain in order to ensure the long-term survival of the human society [78J. 

Even more transparent is the explanation for the existence of negative carrying capacity 
for financial and economic societies, when the variable x represents not population, but 



28 



capitalization. In these cases, negative capacity is nothing but the borrowed resources that 
have to be returned back to the lender. Due to this leverage effect resulting from borrowing, 
a firm can exhibit a fast development. But borrowing cannot last forever. If the firm, 
society, or country does not produce enough and is not able to pay debts, its creditors will 
lose trust and may require early reimbursement, or will simply refuse rolling over the debts, 
as occurred for Greece in May 2010 and Ireland in November 2010. This situation can be 
captured by assuming the existence of a maximum level of debt, beyond which the society 
or country becomes highly unstable due to feedbacks resulting from market forces. Actually, 
Reinhart and Rogoff [79] have recently documented the existence of a strong link between 
levels of debt and countries' economic growth over the last two centuries: Countries with 
a gross public debt exceeding about 90% of annual economic output tended to grow a lot 
more slowly and to exhibit larger default risks. 

Assuming the existence of a maximum debt level beyond which instabilities appear leads 
to the existence of a time t cras h beyond which a crash or at least strong turbulence can occur. 
This is highly reminiscent of the scenario leading to the "great recession" that started in 
2007 worldwide [80J. The minimum crash time £ cms /i is thus given by the condition that the 
debt, represented by the negative carrying capacity, reaches the value 

y(t C rash) = 1 - \b\x(t crash - r) < . (140) 

The crash happens before the critical divergence time t c , 

tcrash *^ j (-^l) 

where the firm or country capitalization is still finite. In such a regime, the accelerated 
growth, fueled by borrowing, leads to a boom that is not supported by increasing productiv- 
ity. This can therefore be called a bubble |28j. As the bubble develops, it eventually reaches 
a threshold level beyond which it becomes unstable, and can therefore be followed by a crash 
at times between t cras h and t c . 

11.3 Boom and crash in society with loss and competition 

The same interpretation as above is applicable for the society with loss and competition, as in 
Sees. 6.4 and 6.6. There, the finite-time singularity arises under a high level of destruction, 
when the destruction coefficient b < — 1 and the initial carrying capacity is negative. The 
hyperbolic divergence occurs at a critical time t c . In Sec. 6.6, the situation is similar to 
that discussed above. The difference between Sec. 6.4 and Sec. 6.6 is that in Sec. 6.4 the 
divergence is defined by the equation 

y(t c ) = 1 - |&|x(* c - r) = . (142) 

Again, a society or a firm with loss and competition, actually, does not reach the point 
of divergence, but becomes bankrupt before this. The fast growth is due to exploiting and 
destroying the carrying capacity. But, after destruction has taken place and reached an 
unbearable level, the boom is followed by a crash. 
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11.4 Species extinction under mutual parasitic symbiosis 



In the parasitic symbiosis of two species considered in Sec. 9.5, there occurs a finite-time 
singularity. Thus, for the symbiotic parameters b < g < 0, the initial carrying capacities are 
such that 

y x (0) = 1 - \b\x z < , y 2 (0) = 1 - \g\x z > . (143) 

For the opposite case, when g < b < 0, the situation is symmetric. Thence, below we shall 
treat the case of Eq. (I143j) without loss of generality. The divergence appears at the critical 
time t c given by the equation 

yi (t c ) = 1 - \b\x(t c )z(t c ) . (144) 

At this time, the population of the species x tends to infinity, while that of the species z 
goes to zero. 

Of course, no realistic population can rise to infinite values. Such a divergence happens 
because of the mutual parasitic symbiotic relations, resulting in the formal appearance of a 
negative effective capacity. As in the cases above, the divergence can be avoided by limiting 
the carrying capacity by a fixed level. This implies that the rise of a parasitic species 
continues only up to some limiting carrying capacity threshold 

y-i i^crash) 1 \b\x{t cras fo)z(t cras i l 

after which the species x dies out by a fast process of extinction at the crash time t cras h < t c . 



11.5 Species extinction under parasitic symbiosis without direct 
interactions 

A finite-time singularity also appears in the case of symbiosis without direct interactions, as 
in Sec. 10.5. This happens when at least one of the species is parasitic. For example, for 
the case b < and b < g. Below, we shall consider this case, since the situation with g < 
and g < b is symmetric. 

When b < 0, this means that the species z is parasitic and destroys the carrying capacity 
of the species x. The finite-time singularity occurs if, at the initial moment of time, the 
effective carrying capacity of species x is negative, 

2/i(0) = 1 - \b\z < , (146) 

while the carrying capacity of species z, 

y 2 (0) = l + ^Xo, (147) 

can be positive or negative, depending on the values of g and x . The divergence of x occurs 
at the critical time t c , where the effective capacity of species x is negative, 

y x {t c ) = 1 - \b\z(t c ) < . (148) 

The population of species z at the moment t c is finite. 
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In the same way as in the previous cases, we understand that this divergence cannot be 
real and there should exist a limiting carrying capacity 



at which the population x is to be set to zero, implying its extinction caused by the parasitic 
species z. This extinction happens at the crash time t CTaa h < t c . 

In all these cases for which there arises a finite-time singularity, it is possible to exclude 
the formal divergence by limiting the carrying capacity to a minimal value y(t cra sh), that is, 
a maximal absolute value \y{t cras h)\. This limiting value can be interpreted as a threshold 
for a change of regime. The overall dynamics, thus, starting with the fast growth of the 
population (or capitalization), is followed by its drop to zero at the crash time t cras h, before 
the critical time t c . 



In this paper, we have suggested a general approach for describing the evolution of popula- 
tions, whose activities influence their carrying capacities. In order to take into account this 
influence, the carrying capacities are to be defined as functions of the society populations. 
This includes the action of a population on its own carrying capacity. In general, the ac- 
tions of populations on the carrying capacities can be delayed, since such actions, generally, 
require time for their realization. 

The approach is illustrated by analyzing the time evolution of a society that acts on 
its own carrying capacity, either by producing the increase of the capacity or by destroying 
it. Different types of societies have been studied, depending on the balance between gain 
and loss and between competition and cooperation. A detailed classification of admissible 
dynamic regimes has been given. 

Two kinds of extreme events have been found to arise, when the society destroys its 
carrying capacity. One is a finite-time death at a death time td and another is a finite- 
time singularity at a critical time t c . The finite-time death describes the extinction of the 
population because of the destruction of the carrying capacity. The finite-time singularity 
signals that the society becomes unstable and its stabilization requires changing the society 
parameters and a transfer to another dynamic regime. The divergence can be avoided by 
limiting the carrying capacity and interpreting the effect as a fast rise of the population 
(or capitalization), followed by its sharp drop. For economic and financial societies, the 
fast growth is understood as a boom or bubble, due to the leverage effect induced by over- 
indebtedness, after which a crash occurs. 

The suggested approach is also illustrated by considering the symbiosis of several species. 
This approach allows us to give a general classification of different symbiosis types. The 
case of two species is analyzed in detail. Extreme events arise when at least one of the 
species is parasitic, destroying the carrying capacity of the other species. Again, there 
can exist two kinds of such extreme events, finite-time death and finite-time singularity. 
Their interpretation is analogous to that given for the case of the self-destructing population 
activity. 

As a general conclusion valid for the different considered situations, we have to say that 
any destructive action of populations, whether on their own carrying capacity or on the 




(149) 



12 Conclusion 
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carrying capacities of co-existing species, can lead to the instability of the society that is 
revealed in the form of the appearance of extreme events, finite-time extinctions or booms 
followed by crashes. 
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Figure 2: Temporal behavior of solutions to equation f )38|) for the initial condition xq = 1, 
lag r = 30, and parameters b = 2 (solid line) and 6 = 3 (dashed-dotted line). Solutions 
x(t) — > oo, when t — )■ oo. 
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Figure 3: Dynamics of solutions to equation ( |38]) for the lag r = 20, symbiotic parameter 
b = 0.5, and the initial conditions xq = 0.1, (solid line) and xo = 1 (dashed-dotted line). 
Solutions x(t) — > x*, monotonically growing by steps, when t — > oo. The stationary point is 
x* = 1/(1 - b) = 2. 



39 



x(t) 




\ . x = 3 

* / 


x = 5 


\ I 


/ 


~V / v...,. 







'o 20 40 60 80 j 100 



Figure 4: Evolution of solutions to equation ( )38|) for the lag r = 20, parameter b = 0.5, and 
initial conditions xq = 3, (solid line) and xq = 5 (dashed-dotted line). Solutions x{t) — > x*, 
monotonically diminishing by steps, when t — > oo. The stationary point is x* = 1/ (1 — 6) = 2, 
the same as in Fig. 3. 



40 



1.25 



0.75 




0.25 



Figure 5: Temporal behavior of solutions to equation f )38|) for the lag r = 20, parameters 
b = —1/3, xo = 1.5 (solid line), and b = —2/3, x = 0.1 (dashed-dotted line). Solutions 
x(t) — > x*, oscillating by steps, when t — > oo. The stationary points are x* = 0.75 and 
x* = 0.6, respectively. 
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Figure 6: Temporal behavior of the solution to equation (j38p for the parameter b = —2 < 
— 1, lag t = 1.18 < t , where r = 1.2092 is defined by fj4"2~|) . with the initial condition 
x = 0.25 < l/|a;o|- Solution x(t) — >■ x* — 1/3, oscillating, when t — > oo. 
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Figure 7: Dynamics of the solution to equation fl38|) for the parameter b = —2 < —1, lag 
Tq < r = 1.5 < 7~i, where r = 1.2092 is defined by f l4"2]) . n = ri(6) w 1.65, with the initial 
condition x = 0.25 < l/\b\. Solution x(t) oscillates without convergence, when t — > oo. 
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Figure 8: Dynamics of the solution to equation (138|) for the parameter b = —2 < — 1, lag 
t = 1.7 > Ti, where T\ ~ 1.65 is defined numerically, with the initial condition x = 0.25 < 
l/|fe|. Solution exists only till the moment t = ~ 19.975, defined by (j49|) . 
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Figure 9: Behavior of the solution to equation f )38|) in logarithmic scale for the parameter 
b = —2 < — 1, lag t = 1.7 > Ti, where T\ ~ 1.65 is defined numerically, and the initial 
condition x — 1 > l/\b\. Solution x(t) diverges, x(t) — > oo, when t — > t c ~ 0.69315. 
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Figure 10: Behavior of the solution, shown in logarithmic scale, to equation ( 1381) for the 
parameter b = —2 < —1, lag r = 0.415 < r c , where r c m 0.416 is defined numerically, with 
the initial condition xq = 1 > l/\b\. Solution x(t) exponentially grows, when t — > oo. 
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Figure 11: Evolution of the solution to equation (155]) for the parameter b = — 1.5 <— 1, lag 
r = 10, with the initial condition x = 1.5 < l/(\b\ — 1). Solution x(t) monotonically decays 
to finite-time death at t d 31.884, defined by (149|) . 
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Figure 12: Stability region for the loss-competition case defined in (J62 



48 



0.8 -I 



0.6- 



0.4- 



0.2- 



I x(t) 




1 * 
1 I 

1 \ 

I \ 
\ 1 


■ 


■ \ * 

\ \ 

\ \ 

\ \ b = - 0.5 
\ \ / 


■ 


\ \ / 
\ v / 

\ / v 


b = 8 


\ \ y 
\t ,/ 





t 



Figure 13: Temporal behavior of the solution to equation (163|) for the lag r = 10, initial 
condition x = 1, with the parameters 6 = —0.5 (solid line) and 6 = 8 (dashed-dotted line). 
The initial condition x < Vl^l f° r ^ < 0. Solutions x(t) monotonically converge to their 
stationary state x* = 0. 
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Figure 14: Solutions x(t) to Eq. (|63|) as functions of time for the parameters r = 0.4 (solid 
line) and r = 0.5 (dashed-dotted line), where r < r = 0.505951, with r defined by (1681) . 
Other parameters are: xq = 1, and b = —2.5. The solutions x(t) converge by oscillating 
towards their stationary point x\ = —1/(1 + b) = 2/3 as t — > oo. 
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Figure 15: Behavior of the solution x(t) in logarithmic scale to Eq. fl63|) as a function of 
time for the parameters b = —2.5, r = 0.6272, and the initial history x = 1- The solution 
x(t) ->• oo as t ->■ t c — 15.1498, where t c is defined by (JTl"]) . 
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Figure 16: Temporal behavior of solutions x(t) to Eq. f )63|) as functions of time for the 
parameters b = —2.5 and r = 0.25 (solid line); r = 0.23 (dashed-dotted line); and the initial 
history xq = 1. The solution x[t) — > x* — 2/3, when t — > oo. 
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Figure 17: Behavior of solutions x(t) to Eq. ( 1791) as functions of time for b = —2.5, lag 
r = 1.5, history Xo = 1 (solid line), and lag r = 0.5, history xq = 0.5 (dashed-dotted line). 
At the moment t = td, defined as 1 + bx(td — r) = 0, solutions monotonically decay to zero, 
x{td) = 0, x{t)\ t =t d = — oo. The death time for the population, represented by the solid line, 
is td = 2.14158, and the death time for the population shown by the dashed-dotted line, is 
t d = 0.58004. 
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Figure 18: Region of stability (shaded) in the parameter plane h — g for the stationary 
solutions in the case of symbiosis with mutual interactions. 
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Figure 19: Comparison of the symbiotic solutions x(t) (solid line) and z(t) (dashed-dotted 
line), for the symbiosis with mutual interactions (j!04p . with the solutions x(t) = z(t) (dashed 
line) of the decoupled equations fillip for the same initial conditions Xq — z — 0.1 < 1, with 
different symbiotic parameters b and g: (a) b = 0.25, g = 0.1 < g c = 0.25, the stationary 
points of the symbiotic equations being x* = 1.411, z* = 1.164; (b) b — 2, g — —0.5, the 
fixed points of the symbiotic equations being x* = 3.562, z* = 0.360; (c) b = —1, g = 2, the 
symbiotic fixed points being x* = 0.293, z* = 2.414; (d) b = —1, g = —2, the symbiotic fixed 
points being x* = 0.707, z* = 0.414. 
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Figure 20: Behavior of the solutions in logarithmic scale to equations (I104p in the case 
of finite-time death and singularity, x(t) (solid line) and z(t) (dashed-dotted line), for the 
parasitic relations with the symbiotic coefficients b = — 1, g = —2, under the initial conditions 
.tq = 1.8, zq = 0.5. For these parameters, the critical time is t c = 0.87. 
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Figure 21: Region of stability (shaded) in the parameter plane b — g for the fixed points in 
the case of symbiosis without direct interactions. 
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Figure 22: Nonmonotonic convergence to stationary states of x(t) (solid line) and z(t) 
(dashed-dotted line) as functions of time, in the case of symbiosis without direct inter- 
actions f 1 1 2 j) : (a) for the initial conditions Xq = 1.8, z = 0.01 and the parameters b = 2, 
g = 0.25. The functions x(t) — > x* and z(t) — > z* , when t — > oo, the fixed points being 
x* = 6, z* = 2.5; (b) for the initial conditions xq = 0.001, zo = 1.5, and the parameters 
b = —0.5, g = 2. Functions x(t) — > x* = 0.25 and z{t) — > z* — 1.5, when t — > oo. 
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Figure 23: Behavior of the solutions x(t) (solid line) and z(t) (dashed-dotted line) in loga- 
rithmic scales in the case of symbiosis without direct interactions ( I120p in the presence of 
the finite-time singularity for the symbiotic coefficients b = —0.75, g = —0.25, and the initial 
conditions xq = 10, zq = 5. The critical time is t c = 0.55303. 
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Figure 24: Finite-time death in the case of symbiosis without direct interactions. Temporal 
behavior of solutions x(t) (solid line) and z{t) (dashed-dotted line) for different symbiotic 
parameters and initial conditions: (a) b = —0.75, g = —0.5, xo = 0.8, zq = 3, the death time 
being t d = 0.204; (b) b = -1.5, g — 1, x — 1, z — 0.1, with the death time t d = 2.412. 
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